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fSJ ' Abstract. The April 1997 giant flare of Mkn 501 is modelled by a stationary 

O^l , Synchrotron-Proton-Blazar model. Our derived model parameters are consistent with 

X-ray-to- TeV-data in the flare state and difl'usive shock acceleration of e~ and p in a 
Kolmogorov/Kraichnan turbulence sprectrum. While the emerging pair-synchrotron 
f ^ ' cascade spectra initiated by photons from 7r°-decay and electrons from tt — > n —* e - 

\Q ', decay turn out to be relatively featureless, /i* and p synchrotron radiation and their 

^^ ' cascade radiation produce a double-humped spectral energy distribution. For the 

^— ^ , present model we find p synchrotron radiation to dominate the TeV emission, while the 

—J. ' contribution from the synchrotron radiation of the pairs, produced by the high energy 

Qs , hump, is only minor. 

O I THE CO- ACCELERATION SCENARIO 

^: 

f^ ■ With its giant outburst in 1997, emitting photons up to 24 TeV and 0.5 MeV in 

K^ ■ the 7-ray and X-ray bands, Mkn 501 has proved to be the most extreme TeV-blazar 

KA \ observed so far (e.g. Catanese et al 1997, Pian et al 1997, Aharonian et al 1999). 

H \ In this paper, we consider the April 1997 flare of Mkn 501 in the light of a 

■ ■ ■ ' modified version of the Synchrotron Proton Blazar model (SPB) (Mannheim 1993), 

and present a preliminary model fit. 

In the model, shock accelerated protons (p) interact in the synchrotron photon 
field generated by the electrons (e~) co-accelerated at the same shock. This scenario 
may put constraints on the maximum achievable particle energies. 

The usual process considered for accelerating charged particles in the plasma jet 
is diffusive shock acceleration (see e.g. Drury 1983, Biermann & Strittmatter 1987). 
If the particle spectra are cut off due to synchrotron losses, the ratio of the maximum 
particle energies 7p,max/7e,max can be derived by equating tacc,p/tacc,e = tsyn,p/tsyn,e, 
with tsyn,p and tsyn,e being the synchrotron loss time scales for p and e~ , respectively. 
We find that for shocks of compression ratio 4 (see Miicke & Protheroe 1999 for a 
detailed derivation) 
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where the "="-sign corresponds to synchrotron loss, and the "<"-sign to adiabatic 
loss determining the maximum energies. 5 is the power law index of the magnetic 
turbulence spectrum {5 = 5/3: Kolmogorov turbulence, 6 = 3/2: Kraichnan tur- 
bulence, and (5 = 1 corresponds to Bohm diffusion). r]e,ma.x is the mean free path at 
maximum energy in units of the particle's gyroradius and F{6, r^cmax) takes account 
of the shock angle 9 (Jokipii 1987). The ratio -F(6', ?7e,max)^e,max/-^(6', r7p,max)^p,max 
can be constrained by the variability time scale tvar, requiring t^arD > tacc,p,max {D 
= Doppler factor, tacc,p,max = acceleration time scale at maximum particle energy) 
for a given parameter combination. As an example, we adopt D = 10, B = 20G 
and tvar = 2 days. Eq. 1 then restricts for these parameters the ratio of the al- 
lowed maximum particle energies to the range below the solid lines shown in Fig. 1. 
Points exactly on this line correspond to synchrotron-loss limited particle spectra 
which are accelerated with exactly the variability time scale. 
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FIGURE 1. Allowed parameter 

space (shaded area) for 7p,max: 
7e.max in the SBP-model for typi- 
cal TeV-blazar parameters (B=20 G, 
D=10, ui = 0.5c, /3 == 1, 
tvar — 2days) and for different mag- 
netic turbulence spetra I{k) oc k^ . 
The diamond symbol corresponds to 
the Mkn 501-model presented below. 
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In hadronic models vr photoproduction is essential for 7-ray production. The 



threshold of this process is given by emaxTp.max = 0.0745 GeV where ej, 
the maximum photon energy of the synchrotron target field. Inserting e^ 
3/87emax-^/(4-414 X lO^^G) mgC^ into the threshold condition, we find 
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which is shown in Fig. 1 as dashed line for various magnetic field strengths. Together 
with Eq. 1, the allowed range of maximum particle energies is then restricted to 
the shaded area in Fig. 1. 



II THE MKN 501 FLARE IN THE SYNCHROTRON 
PROTON BLAZAR (SPB) MODEL 

We assume the parameters used in Fig. 1, and that the co- accelerated e~ produce 
the observed synchrotron spectrum, unhke in previous SPB models, and this is the 
target radiation field for the j97-interactions. This synchrotron spectrum, and its 
hardening with rising fiux, has recently been convincingly reproduced by a shock 
model with escape and synchrotron losses (Kirk et al 1998). We use the Monte- 
Carlo technique for particle production/cascade development, which allows us to 
use exact cross sections. 

FIGURE 2. Mean energy loss 
time of p for synchrotron radiation 
(syn), TT-photoproduction (tt) and 
Bethe-Heitler pair production (e ), 
and for tt - and /x -synchrotron ra- 
diation (syn TT, syn /i) for B — 19.6 G 
with their mean decay time scales (de- 
cay TT, decay /z) in the SPB-model. 
The acceleration time scale (ace), 
based on Kolniogorov turbulence, is 
calculated for ui — 0.5c, 7]p = 40 and 
shock angle 9 = 85°. Its curvature 
reflects the influence of the shock an- 
gle. The adiabatic loss time (adiab) 
is assumed to be R/ui « Dtyar- AH 
quantities are in the jet frame. 
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For simplicity we represent the observed synchrotron spectrum (target photon 
field for the p7-collisions) as a broken power law in the jet frame with photon power 
law index 1.4 below the break energy of 0.2 keV, and index 1.8 up to 50 keV. 

The variability time scale restricts the radius R of the emission region. For our 
model we use t^ar.x ~ 2 days (Catanese et al 1997), and find R ~ 2.6 x lO^^cm for 
D = 10, B = 19.6 G. With these parameters the 77-pair production optical depth 
reaches unity for ^ 25 TeV photons. 

Our model considers photomeson production (simulated using SOPHIA, Miicke 
et al 1999), Bethe-Heitler pair production (simulated using the code of Protheroe & 
Johnson 1996), p synchrotron radiation and adiabatic losses due to jet expansion. 
The mean energy loss and acceleration time scales are presented in Fig. 2. 

Synchrotron losses, which turn out to be at least as important as losses due to 
71 photoproduction for the assumed 2 day variability, limit the injected p spectrum 
oc 7~^ to 2 < 7p < 4.4 X 10"*^°. This leads to a p energy density Up ~ 0.2 TeV/cm , 



which is bracketed by the photon energy density Ut, 



argct 



0.01 TeV/cm , and a 



magnetic field energy density ub ~ 9.5 TeV/cm . With ub ^ ^target significant 
Inverse Compton radiation from tlie co-accelerated e~ is not expected. 

Raclien & Meszaros (1998) noted the importance of synchrotron losses of fi"^- 
(and TT^-) prior to their decay in AGN jets and GRBs. For the present model, the 
critical Lorentz factors 7^^ ~ 3 x 10^ and 7^ ?a 4 x 10^°, above which synchrotron 
losses dominate above decay, lie well below the maximum particle energy for /x^, 
while vr^-synchrotron losses can be neglected due to the shorter decay time. 
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FIGURE 3. Left: Average emerging cascade spectra initiated by 7r°-decay (solid line) and 
TT -decay synchrotron photons (dashed-dottcd line). Right: Average emerging cascade spectra 
initiated by p- (solid line) and y^^-synchrotron photons (dashed-dotted line). 



The matrix method (e.g. Protheroe & Johnson 1996) is used to follow the pair- 
synchrotron cascade in the ambient synchrotron radiation field and magnetic field, 
developing as a result of photon-photon pair production. The cascade can be initi- 
ated by photons from 7r°-decay ("vr^-cascade"), electrons from the n^ -^ fi^ -^ e^- 
decay ("Tr'^-cascade"), e^ from the proton-photon Bethe-Heitler pair production 
( "Bethe-Heitler-cascade" ) and p and yU-synchrotron photons ( "p-synchrotron cas- 
cade" and "/i^-synchrotron cascade"). In this model, the cascades develop linearly. 

Fig. 3 shows an example of cascade spectra initiated by photons of different ori- 
gin, and for the parameter combination given above. 7r°- and vr^-cascades obviously 
produce featureless spectra whereas p- and yU^-synchrotron cascades cause the typ- 
ical double hump shaped SED as observed in 7-ray blazars (see also Rachen, these 
proceedings). The contribution from Bethe-Heitler cascades turns out to be neg- 
ligible. Direct p- and /x^-synchrotron radiation is responsible for the high energy 
peak, whereas the low energy hump may be either synchrotron radiation from the 
directly accelerated e^ and/or by pairs produced by the "low energy hump" . 

Adding the four components of the cascade spectrum in Fig. 3 and normalizing 
to an ambient, accelerated p density of ntot,p = 7 cm~^, we obtain the SED shown 



in Fig. 4 where it is compared with the multifrequency observations of the 16 April 
1997 flare of Mkn 501. 
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FIGURE 4. Present model (histogram) in comparison with the data of the 16 April 1997-flare 
of Mkn 501. Photon absorption on the cosmic diffuse background radiation field is not included 
in the model. Straight solid lines: parametrization of the observed, curved synchrotron spectrum 
(BeppoSAX & OSSE) by Bednarek & Protheroe (1999) and observed TeV-emission corrected for 
cosmic background absorption (Bednarek & Protheroe 1999) with peak energy output ^ 2 TeV; 
the 100 MeV upper limit is from Catanese at al 1997 (observed 9-15 April 1997), diamonds: nearly 
simultaneous (uncorrected) Whipple data (Catanese et al 1997); dashed-dotted line: synchrotron 
target spectrum. 
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